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Abstract 

In this paper, we develop an analytical approach to Doppler cooling of atoms by one- or two- 
photon transitions when the natural width of the excited level is so small that the process leads to 
a Doppler temperature comparable to the recoil temperature. A "quenching" of the sharp line is 
introduced in order to allow control of the time scale of the problem. In such limit, the usual Fokker- 
Planck equation does not correctly describe the cooling process. We propose a generalization of the 
Fokker-Planck equation and derive a new model which is able to reproduce correctly the numerical 
results, up to the recoil limit. Two cases of practical interest, one-photon Doppler cooling of 
strontium and two-photon Doppler cooling of hydrogen are considered. 
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I. INTRODUCTION 



Laser cooling is one of the greatest achievements of Atomic Physics in the few last decades, 
having open a thorough new field of applications and led to the first experimental demon- 
stration of the Bose-Einstein condensation in dilute gases |^. The simplest laser cooling 
process is the well-known Doppler cooling by a red-detuned standing wave interacting with 
a two-level atom . A usual way of describing Doppler cooling is to derive a Fokker-Planck 
equation, leading to a limit "Doppler temperature" given by To = hV /{2kB), where V is the 
width of the transition and fc^ is the Boltzmann's constant. The Fokker-Planck equation is 
valid as long as the mean velocity of the atoms is large compared to the atomic velocity shift 
due to a fluorescence cycle, which is of the order of the recoil velocity Vr = hk^/M, where 
ko = uo/c is the wavenumber of the atomic transition and M is the mass of the atom. The 
Doppler temperature should thus be much larger than the recoil temperature T^. = Mv'^/ks, 
which implies 

r > koVr . (1) 

(Note that k^Vr is the Doppler effect associated to the recoil velocity). In the present paper, 
we shall consider two cases in which the inequality ([l|) is not satisfled, namely the one-photon 
Doppler cooling of strontium and the two-photon Doppler cooling of hydrogen. 

Doppler cooling of strontium has recently been studied by experimental groups, with aims 
varying from achieving high phase-space densities [Q] to the study of the coherent backward 
scattering of light by a disordered cold-atom cloud The idea is to take advantage of 
the very sharp (F ^ 7.5 kHz), spin-forbidden, line ^5*0 — ^ Pi whose theoretical Doppler 
temperature is comparable to the recoil limit. 

Laser cooling of hydrogen atoms, on the other hand, is still a goal to be achieved. Laser 
cooling of trapped hydrogen atoms was experimentally demonstrated , but cooling of free 
atomic hydrogen is still a defy for experimentalists. The main difficulties are that i) hydrogen 
forms tightly bonded molecules that should be dissociated before any experiment on atoms 
can be done, and ii) the optical transitions starting from the fundamental state all fall in the 
vacuum-ultraviolet (VUV), e.g. at 121 nm for the 15* — 2P transition. Moreover, although 
Bose-Einstein condensation of hydrogen was achieved without use of such techniques, 
laser cooling of hydrogen would certainly improve and simplify the procedure leading to 
condensation, which can be of considerable practical relevance. Cooling of atomic hydrogen 
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by pulsed two-photon transitions was suggested in 1993 by M. Allegrini and E. Arimondo 
using TT pulses on the 15* — 3S transition (wavelength of 200 nm for two-photon transitions). 
Two-photon transitions present the advantage of dividing the laser frequency by two, which 
means that the required radiation falls in more favorable wavelength domains. The recent 
progress on the generation of CW, VUV, laser radiation, in the context of metrological 
studies of the IS — 2S transition allows one to realistically envisage the two-photon 
Doppler cooling of hydrogen in the continuous wave regime on the very sharp IS — 2S 
transition (the lifetime of the 25" state is 0.2 s), as we suggested in a previous work @|. 

Laser cooling relies on the ability of atoms to perform a great number of fluorescence 
cycles in which momentum is exchanged with the radiation field in a short time (compared 
e.g. to the collision time). From this point of view, sharp lines are not suitable for cooling, 
because the cooling time becomes very long. On the other hand, the Doppler temperature 
is proportional to the linewidth of the excited level (this result is also valid for two-photon 
transitions); from this point of view, sharp lines are interesting for cooling. In order to 
conciliate these antagonistic characteristics of sharp transitions, we proposed in Ref. 
to use the "quenching" of the 25" state of hydrogen to control the cycling frequency of the 
process and we shown that this process allowed to approach the recoil temperature where 
the Fokker-Planck equation breaks down. In the present work we develop a refined analysis 
of the Doppler process valid up to the recoil limit both for one- and two- photon Doppler 
cooling. 



II. DOPPLER COOLING WITH ONE- AND TWO-PHOTON TRANSITIONS 

Let us consider an atom of mass M and velocity v parallel to the z-axis (Fig. |1|) interacting 
with two counterpropagating waves of angular frequency ujl with puJi = uq + S, {p = 1,2 for 
resp. one- and two-photon transitions), ujq/2tt = 2.5 x 10^^ Hz for H and 4.4 x 10^*^ for Sr). 
We have also = Uq/c = pk^. We consider three atomic levels, namely the ground state 
\g)^ which is coupled by the radiation, via one- or two-photon transitions, to a sharp level 
\s) of linewidth F^, and the wide level \w) which is coupled to the sharp level |s) through 
a quenching process of rate Vq [10| and to the ground- level \g) by spontaneous emission 
(natural linewidth F^;). We introduce two simplifying assumptions: i) The wide level \w) 
can be eliminated adiabatically, so that everything happens as if the sharp level has, due to 



3 




FIG. 1: Atomic levels involved in the cooling process. |s) is the sharp level and \w) is the wide 
level. These levels are coupled by a quenching process of rate Tq. Spontaneous emission from the 
level \s) is very small. 

the quenching process, an effective hnewidth 

r' = r.TT^ = ■ (2) 

where is a controllable quenching ratio, ii) We consider that the sharp and wide levels 
correspond to the same transition frequency with respect to the ground state. This is a 
very good approximation for the case of hydrogen, where the quenching couples levels 25" 
and 2P, that are quasi-degenerated. In the case of the strontium the transition frequencies 
differ by about 30 %, but this does not affect qualitatively our conclusions. The velocity 
shift corresponding to spontaneous emission of a photon is thus = hko/M 3.1 m/s, or 
~ 1-2 mK for H, and I'r ~ 7 mm/s and ^ 0.53 /jK, for Sr. 

We shall now introduce some specifics of the two-photon Doppler cooling. Four two- 
photon absorption process are allowed: i) absorption of two photons from the +z- 
propagating wave (named wave "+" in what follows), with a rate Fi and corresponding 
to the a total atomic velocity shift of -1-^^; ii) absorption of two photons from the —z- 
propagating wave (wave "— "), with a rate F_i and atomic velocity shift of —Vr', Hi) the 
absorption of a photon in the wave followed by the absorption of a photon in the wave 
"— " , with no velocity shift and iv) the absorption of a photon in the wave "— " followed by 



the absorption of a photon in the wave "+", with no velocity shift. The two latter process 
are indistinguishable, and the only relevant transition rate is that obtained by squaring the 
sum of the amplitudes of these process (called Tq). Also, these process are "Doppler-free" 
(DF) as they are insensitive to the atomic velocity (to the first order in v/c) and do not shift 
the atomic velocity. Thus, they cannot contribute to the cooling process. As atoms excited 
by the DF process must spontaneously decay to the ground state, this process heats the 
atoms. In the limit of low velocities, the transition amplitude for each of the four processes 
is the same. One thus expects the DF transitions to increase the equilibrium temperature 
by a factor of two, which is verified both numerically and analytically [^]. Note that the 
analysis of the one-photon process can be deduced straightforwardly from the results below 
by setting the Doppler-free term to zero. The transition rates are given by 



r» = 1^ (3b) 

where p = 1,2 corresponds to, resp., one- and two-photon transitions. T±i describe, respec- 
tively, the absorption from the and "— " waves, and Fq the DF transition rate (note that 
this terms is equal to zero for one-photon transitions). / = I /Is where Ig is the saturation 
intensity, 6 is the detuning divided by F^, K = k^Vr/T^ ^ 0.26 for H and K ^ 3.2 x 10~^ 
for Sr, and V = v/vj.. The only qualitative difference between one- and two-photon Doppler 
cooling, in the present context, is the presence of Doppler-free terms in the latter case. Note 
that, with the above normalizations, the usual Fokker-Planck approach is expected to break 
down when g < \6\ K. 

III. DERIVATION OF THE GENERALIZED FOKKER-PLANCK EQUATION 

The rate equations describing the evolution of the velocity distribution n(y, t) and n*(y, t) 
for, respectively, atoms in the ground and in the excited level are 

^^^^ = - [F_i(V) + Fo + T,{V)\n{V, t) + [n\V -l)+n*{V+ 1)] (4a) 
Bn*(V t) 

= F_i(V^ - l)n{V - 1, t) + Ton{V, t) +T,{V + l)n{V + 1, t) - rn*{V, t) . (4b) 

ot 
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The deduction of the above equations is quite straightforward (cf. Fig |l]). The first term 



in the right-hand side of Eq. ( pi]) describes the de-population of the ground-state velocity 
class V by two-photon transitions, whereas the second term describes the re-population of 
the same velocity class by spontaneous decay from the excited level. In the same way, the 
three first terms in the right-hand side of Eq. ( ^ ) describe the re-population of the excited 
state velocity class V by two-photon transition, and the last term the de-population of this 
velocity class by spontaneous transitions. For each term, we took into account the velocity 
shift {V V ±1) associated with each transition and supposed that spontaneous emission 
is symmetric under spatial inversion. 

For laser intensities below the saturation intensity of the transitions, one can adiabatically 
eliminate the population of excited level and reduce the Eqs. (|^) and(|4ED to one equation 
describing the evolution of the ground-state population: 



dn{V,t) 
Jt 



H 7: 1 7: — 



n{V,t) + 

^{Fo [n{V -l,t)+ n{V + 1, t)] + F_i(V - 2)n{V - 2, t) 
+ Ti{V + 2)n{V + 2,t)^ (5) 

The set of ordinary differential equations in Eqs. (^al) and (|bD or Eq. (1) has no known 
analytical solution and must be solved numerically. However, it is possible to develop some 
analytical approaches following the general idea leading to the Fokker-Planck equation (FPE) 



for the time-evolution of the velocity distribution (see Refs. 0, |T2[)- Taking the continuous 
limit of Eq. (^ with respect to V, we obtain the following "generalized Fokker-Planck 
equation" (GFPE) 0: 

dn d{fin) d'^ (12^) 



dt dV dV^ 
where the functions fi{V) and f2{V) are given by: 



(6) 



MV) = T^,iV)-T,iV) (7a) 
MV) = T_,{V) + UV) + ^ . (7b) 

In this form, we keep explicitly the ^-dependent coefficients fi{V) and f2{V). Note however, 
that this limit implies a smooth variation of the velocity distribution n{V) as well as of the 
rates r-|-i(V") [on a l^-interval which is 0(1)]. However, the rates r-i-i(l^) present sharp peaks 
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FIG. 2: Velocity distribution for two-photon Doppler cooling obtained from the GFPE [curve (a)], 
and FPE [curve (b)] compared to the numerical simulations of Eq. (||) (dashed line). For the sake 
of comparison all distributions are set to n(V) = 1 for V = 0. The distribution temperatures are 
T = 2.19 for (a), T' = 3.88 for (b) and 2.16 for the numerical simulation (temperature in recoil 
units). Parameters are g = 0.2 and K = 0.26. 



at the values V = ±6/K (in the hmit g ^ 0) and the GFPE is thus not vahd around these 
values. Note that the standard FPE can be obtained from the above equations in the limit 
V -> 0: 



fiiV) 
f2{V) 



2V 



dV 



2T_,iV = 0) + ^ 



(8) 
(9) 



where the derivative is evaluated at = 0. 

To our knowledge, the above system has no analytical solution. However, we can focus on 
their asymptotic steady-state solution. By setting dn/dt = in Eq. (|j), we straightforwardly 
obtain: 



n{V) 



1 



hiv) 



exp 



h{v) 

f2{V) 



dV 



(10) 



The steady-state distribution is then obtained by integration of Eq. ( pUj ) (/1//2 is a ratio of 
polynomials in V and the integration is obtained by standard methods). One gets: 



n{V) 



f2iV) 



exp 



2 (52 + g^) S 



K^{76^ + 2252^2 _ ^4) ^(754 + 22^2^2 _ g4 



tan 
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[provided {7S^ + 22S'^g'^/4 — (7^/16) > or more simply \S\ > g/10. Note that this inequahty 
is fulfilled in most practical situations, since \6\ must be high enough to allow the system 
to distinguish the "+" and the "-" cooling beams]. The solution of the FPE, on the other 
hand, is the gaussian velocity distribution: 

1 \6\K 



n'{V) = exp 



(12) 



2 52 + ^2/4 

These velocity distributions are shown in Fig. § and compared to the results of a direct 
numerical simulations of Eq. (|^). Eq. (pUf ) compares very well to the numerical solution for 
\V\ < \S\/K ~ 4 (that is, for the central peak of cold atoms), while the FPE distribution 
Eq. (|l^) leads to a broader (hotter) distribution. Both approaches fail to correctly describe 
the (uninteresting) background of hot atoms. 

We can deduce the temperature of the distribution from its the central peak of cold atoms 
(which typically corresponds to \V\ < \S\/K, that is, to atoms strongly interacting with the 
radiation) and we neglect the wide background of hot atoms. In this limit, the velocity 
distribution fits to a gaussian shape exp{—V'^/2T) where T is the temperature in units. 
For the GPFE, we obtain the following analytical expression: 

T = (^^ + ^V4)2 

whereas for the FPE distribution the same approach gives: 

52 + o2 /4 

r = ^ ' (14) 

K\5\ ^ ' 

Fig. ^ compares both temperatures to the results of numerical simulation. One observes 
that the values of T' increasingly disagrees with the numerical results as one approaches 
the recoil limit (as expected), and predicts values that are typically too hot by a factor ^ 3 
for detunings close the Doppler minimum \5\ ~ g/2 [obtained from Eq. ([l^)]- For higher 
detunings (|5| ^ K), both the FPE and GFPE approach the same limit T ^ \5\/K. Note 
that for I (5 1 < g/2 the GPFE approach also breaks down, leading to sub-recoil temperatures. 
In this case, the peak of cold atoms disappears into the background of hot atoms, it is then 
hard define a temperature. However, this regime is not of great experimental interest. 
In the case of the cooling by one-photon transitions, we obtain the velocity distribution 
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detuning (T^) 



FIG. 3: Temperature vs. detuning for two-photon cooling. The curve (a) corresponds to the GFPE 
and the curve (b) to the FPE (both represented by sohd hnes). The numerical values corresponds 
to the dashed-dotted line). The agreement is clearly better for the GFPE. Parameters are g = 0.2 
and K = 0.26. 

simply by setting Fq = in Eqs. ( [7b| ) and (H): 



for the GFPE and FPE cases respectively. One deduces the temperatures of the cold central 
peak in the same way as above. The result is exactly T/2 and T'/2. The factor two is a 
consequence of the Doppler-free transitions that increases the equilibrium temperature, as 
we pointed out in Sec. I and discussed in more detail in Ref. 0. 

Figure § shows the distributions compared to the results of a numerical integration of 
Eq. The GFPE is seen to correctly describe the velocity distribution and the temperature 
up to the recoil limit. Note that 6 and g must be of the order of K in order to obtain 
temperatures of the order of the recoil temperature [see Eqs. ([131) and (|14[)]. 

In the case of hydrogen, a stray effect can show up for the high laser intensities: the 
photoionization of the excited state, causing a decreasing in the number of cooled atoms. 
We evaluated the impact of this effect in Ref. 0, and showed that it does not limit the 



n{V) 



/2(^)|ro=o 




(15) 



method. 
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FIG. 4: Velocity distribution for one-photon transitions : GFPE [curve (a)], and FPE [curve 
(b)] (solid lines) and numerical simulation based on Eq. @ (dashed line). The temperatures 
are T = 1.30 for (a), T' = 2.12 for (b) and 1.4 for the numerical simulation. Parameters are 
K = 3x 10"^ for strontium, 5 = 4:K and g = 0.8K. 

IV. CONCLUSION 

In conclusion, we have developed a generalized Fokker-Planck approach allowing to de- 
scribe Doppler cooling on very sharp lines up to the recoil limit, and showed that the resulting 
expressions compare very well to the "exact" numerical. This approach can in principle be 
extended to a more precise description of particular cases. 
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